Geometry-induced pulse instability in microdesigned catalysts: 
the effect of boundary curvature 
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We explore the effect of boundary curvature on the instability of reactive pulses in the catalytic 
oxidation of CO on microdesigned Pt catalysts. Using ring-shaped domains of various radii, we 
find that the pulses disappear (decollate from the inert boundary) at a turning point bifurcation, 
and trace this boundary in both physical and geometrical parameter space. These computations 
corroborate experimental observations of pulse decollation. 

PACS numbers: 82.40.Ck, 82.40.Np, 05.45.-a 



I. INTRODUCTION 



Pattern formation in spatially uniform media has been 
studied extensively (see e.g. the review by Cross and Ho- 
henberg |l| and, specifically for catalytic surfaces, the re- 
cent review by Imbihl 2]). Since the pioneering Photoe- 
mission Electron Microscopy (PEEM) work of Ertl and 
coworkers in 1989 3], high vacuum CO oxidation on sin- 
gle crystal Pt catalysts has been a paradigm for such 
studies. Current research increasingly focuses on the in- 
terplay between spontaneous pattern formation and spa- 
tial medium nonuniformity, as well as with spatiotem- 
poral variations of the medium properties. Intention- 
ally designing the geometry of the medium at the mi- 
croscopic level has been finding an increasing number of 
applications across disciplines in recent years: from the 
guidance of cell migration through micropatterning Q 
and the template-based self-assembly of diblock copoly- 
mers (e.g 5]) to the observation of binary fluids p and 
phase separation in confined geometries (e.g. 7]). Ex- 
amples of complex, microdesigned geometries used in the 
study of reacting systems include Belousov-Zhabotinsky 
(BZ) catalyst microprinting 8] and CO oxidation on mi- 
crolithographically shaped domains on Pt catalysts Q. 
The study of reactive pattern formation in random het- 
erogeneous media is a topic of active current modelling 
work (e.g. OHEl); natural heterogeneous media (such 
as water-in-oil microemulsions in [lg) are also being 
studied. Microfluidics is another field that hinges on the 
micron-scale design of geometry to control flow and trans- 
port (see e.g. the review |l3j): beyond passive geometry 
design, extensive developments in active addressing, for 
example in control of droplet breakup 14] or micromix- 
ing [T3 have been realized. Actively and spatiotempo- 
rally altering the activity of chemically reactive media 
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is exemplified in 16] for thermally addressed heteroge- 
neous catalytic reactions and in [l7l for photosensitive 
BZ media. The recent work of marries microfluidic 
technology with the active addressing of pattern form- 
ing reactive media in a developmental biology context, 
using temperature fields to perturb pattern formation in 
drosophila embryos. 

The interplay of nonlinear reaction-diffusion dynam- 
ics with the effects of boundaries has been extensively 
studied both in theory and in experiments; a remark- 
able early observation was that rotating spiral waves of 
a BZ reagent could be initiated close to domain bound- 
aries with sharp corners [HI- Observations in catalytic 
reactions include that (a) below a certain critical domain 
size, the frequency of rotating waves in a small circular 
domain is affected by the domain size 20]; and (b) that 
the interaction of patterns with inert or active boundaries 
causes the pinning, transmission and boundary backup of 
spirals, and instability of pulses turning around corners 
|2lj . The presence of boundaries can also give rise to 
new patterns that have not been observed in homoge- 
neous reaction-diffusion systems 

Travelling pulses in excitable media constitute an im- 
portant building block of pattern formation. They have 
been documented in a variety of nonlinear dynamical 
systems, including CO oxidation on Pt(llO), arising 
from a localized, finite perturbation of a linearly sta- 
ble steady state [22, Hj| • The instabilities and bifurca- 
tions of solitary travelling pulses in the CO oxidation on 
Pt(llO) have been well studied in the one-dimensional 
case for a wide range of parameters [2J, [23 . A recent 
study explored the effects of anisotropic diffusion on two- 
dimensional pulse propagation in the thin ring limit |2£j . 
However, the experimentally observed pulses in the cat- 
alytic CO oxidation on Pt(llO) propagate on a truly two- 
dimensional surface and interact with boundaries and 
heterogeneities of finite size; a more detailed study of the 
boundary effects on two-dimensional pulse propagation 
has been missing. The use of a simplified two-variable 
reaction-diffusion model [2^ in previous studies of pulses 
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in CO oxidation also creates some difficulty in compar- 
ing computational results with experimental data. De- 
rived from the full, three- variable model [2^ through ap- 
proximations and simplifications, the parameters in the 
two- variable model do not have simple, explicit physical 
meaning; variation of such model parameters is not easy 
to realize experimentally. 

In this paper, we use the full, three-variable Krischer- 
Eiswirth-Ertl (KEE) model for CO oxidation on Pt(llO) 
to study and analyze the behavior of travelling pulses 
in fully two-dimensional ring structures. In particular, 
we perform stability and continuation computations with 
respect to both physical parameters (such as the partial 
pressure of CO (Pco)) and geometric parameters such as 
the ring radii. The curvature of an inert boundary on 
which a pulse is attached plays a crucial role in affecting 
the pulse dynamics; the ring geometry provides an ideal 
setting in which to study this fundamental interaction 
since the (constant) curvature is prescribed and can be 
varied systematically when constructing the domain. 

The paper is organized as follows: We begin with a 
brief introduction to the KEE reaction-diffusion model 
for catalytic CO oxidation on Pt(llO) in Section |H] Sec- 
tion IIIII shows the effect of boundary curvature on the 
stability of two-dimensional travelling pulses in two dis- 
tinct cases: for fixed geometry, the dependence of pulse 
stability on the partial pressure of CO is first studied ; 
then, keeping Pco constant, we study the dependence of 
pulse stability on the boundary curvature — both inner 
and outer — . Representative experimental results sup- 
porting the calculation are presented in Section HVI We 
briefly summarize our results and conclude in Section IVl 



II. MODELING 

We use the KEE reaction-diffusion model for CO oxi- 
dation on Pt(llO) with a surface phase transition [28^ . 
This surface reaction follows a Langmuir-Hinshelwood 
mechanism: 

CO + * ^ CO ads 
2 * +0 2 -> 20 ads 

CO ads + O ads 2 * +C0 2 T 

accompanied by a 1 x 2 ^ 1 x 1 phase transition of the 
Pt(llO) surface due to CO adsorption. The equations in 
this kinetic model are 



w k 3 (f(u) - w) 



(3) 



where v and w denote the surface coverage of CO and 
O, and the surface fraction of the lxl phase respec- 
tively. The adsorption rate constant for CO and O2, k u 
and k v respectively, are set to be constant within the 
temperature range considered in this paper. The rate 
constants /ci,/c2 and ks for the desorption, reaction and 
surface phase transition are given by the Arrhenius for- 
mula ki = (k°)iexp(—Ei/RT); T is the temperature of 
the single crystal. The function f(u) has been fitted to 
experimental data to give the rate of surface phase tran- 
sition as a function of u, the coverage of CO, as follows: 





-1.05n 2 +Q.3n-Q.Q26 
-0.0135 
1 



for u < 0.2 
for 0.2 < u < 0.5 
for u ^ 0.5 



The variable v can, in principle, be adiabatically elimi- 
nated to give an activator-inhibitor type PDE system of 
two equations 

We discretize the three model equations over a two- 
dimensional ring-shaped domain using the Finite Ele- 
ment Method as implemented in the package FEMLAB 
( http : / / www . com sol . conf . FEMLAB provides a dy- 
namic simulator for the problem; we want to find stable 
and unstable stationary (in a travelling frame) states, 
perform continuation/bifurcation calculations, and com- 
pute the linearized stability of the solutions we find. We 
used matrix-free linear algebra algorithms to perform 
these computations, with the right-hand-side of our equa- 
tions provided on demand by FEMLAB. Newton-Krylov 
GMRES (see e.g. [29|) was used for stationary state 
computations, coupled with pseudo-arclength continua- 
tion to trace the bifurcation diagrams; the Jacobian of 
the discretized equations was exported from FEMLAB 
into a MATLAB code for this purpose. ARPACK was 
used within MATLAB for stability computations of the 
resulting generalized eigenvalue problem. 

Due to the symmetry of the domain, we rewrite the 
system of differential equations in a frame which rotates 
at an angular speed uj. The model equations written in 
co-moving frame and polar coordinates are given by: 



u = k u s u pco I 1 — ( — ) ) -kiU-k 2 uv+D u X7 u+udou 



(4) 



u = k u s u pco I 1 - ( — ) ] -k 1 u-k 2 uv + D u V u (1) 
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v = k v po 2 (ws Vl + (l-w)s V2 ) ( 1 
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w = ks(f(u) — w) + ujOqw 



-k2uv j ruodov 
(5) 

(6) 
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The diffusion of CO is chosen to be isotropic for sim- 
plicity. By appending an extra pinning condition to re- 
move rotational invar iance, we can find the appropriate 
co-rotating speed uo as well as the steady pulse shape 
travelling with this speed, as steady state solutions to 
Eqn. 0] and El It is worth noting that ways to find 
the "right" co-rotating speed, even during transient evo- 
lution, have been devised using a "template based" ap- 
proach (IMUJEI). 

We want to study the influence, on the stability of the 
travelling pulses, of the boundary curvature by varying 
the ring radii. The ring radii do not appear explicitly in 
the equations, however; they only appear in the boundary 
conditions. We therefore first carry out a linear transfor- 
mation of coordinates, mapping the physical domain (the 
strip between the inner radius r\ and outer radius r ) to a 
constant computational domain (between the fixed values 
PuPo) 

Let r, 9 be the independent variables in the original 
polar coordinates. The diffusion term in Eqn. is 
expressed as 



1 du 

r dr 



cPu 
dr 2 



1 d 2 u 



, r e [ri,r ], e (0,2tt] 



Define another set of polar coordinates (p, < 



( \ Po Pi . 

(r - n) x h pi 



P £ [pi,Po] 

<j> E (0, 2tt] 



and express the term of CO diffusion in the new polar 
coordinates: 



V 2 u 



yjc\ci du d 2 u C2 d 2 u 
p dp Cl dp 2 p 2 d(j) 2 



ci 
c 2 



Vu 



y/CiC 2 - Ci) 



Vu 
(7) 



where 



ci 



Po ~ Pi 



' O 1 I 



and Vu is the gradient of u in the new polar coordinates. 
Now r Q and T{ appear explicitly in ci, C2- By substituting 
Eqn. in Eqn. QJ, we obtain a new PDE system that 
can be discretized and rewritten in the co-moving frame 
to solve for the steady pulse solutions. 



III. COMPUTATIONAL RESULTS 

A. Influence of the presence of boundary 

We recall some results for simple, one-dimensional soli- 
tary pulses with respect to the variation of Pco, as shown 
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FIG. 1: Bifurcation diagram of ID pulses with respect to the 
partial pressure of CO. The variable c is the propagating speed 
of the steady pulse (a representative one is shown in the inset). 
The turning point (saddle-node bifurcation) is indicated by 
a red circle. Solid (dashed) lines denote stable (unstable) 
steady pulse solutions. T — 540 K, Po 2 — 1-33 x 10~ 4 mbar. 
The inset shows the typical shape of u, v and w at the steady 
state marked by the arrow. 



in Fig. The travelling speed of a stable solitary pulse 
decreases as Pco increases. When Pco is above some 
critical value (4.599 x 10~ 5 mbar in Fig. [I}, the medium 
can not support pulses any more. Checking the leading 
eigenvalues of the linearization of the system equations 
we identify the bifurcation at which the pulse is lost as a 
turning point (or saddle-node) instability: a single eigen- 
value crosses zero and becomes positive as the branch 
of travelling pulses turns around and becomes unstable 
(saddle- type). In a two-dimensional ring-shaped domain, 
the shape and stability of travelling pulses is affected by 
the presence of inert boundaries, as shown in Fig. [21 
The insets show the oxygen coverage profile correspond- 
ing to the steady pulse at different locations along the 
diagram. The inner and outer radius of the ring are fixed 
at dimensionless values of 20 and 30 respectively (1 unit 
corresponds to approximately 3.8 pms). Comparing the 
bifurcation diagrams in Fig. [21 and [I] shows that the turn- 
ing point bifurcation is still present, but now we find that 
the pulse loses stability before the turning point due to a 
Hopf bifurcation; simulations indicate that this bifurca- 
tion is subcritical. We plot the eigenmodes correspond- 
ing to the leading eigenvalues at points H and p along 
the diagram in Fig. [3] The leading eigenmodes at the 
Hopf bifurcation point H, Fig. EJb) and|3fc), constitute 
a complex pair crossing the imaginary axis and render- 
ing the pulse unstable. The spatial structure of these 
eigenmodes clearly resides close to the inner boundary, 
suggesting that this boundary is responsible for the in- 
stability. An eigenvalue zero corresponding to the rota- 
tional invariance is also found; its eigenvector (Fig. EJa) 
and|3fd)) can be obtained from the pulse shape by differ- 
entiation with respect to 6. When the pulse propagates 
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FIG. 2: Bifurcation diagram of pulses in a 2D ring showing 
the influence of the CO partial pressure on the angular speed 
(ui) of the travelling pulse at stationarity. Point if, marked 
with a red circle, denotes a Hopf bifurcation. Solid (dashed) 
lines denote stable (unstable) steady pulse solutions. The 
insets show profiles of O coverage. The color is scaled from 
blue to red according to the coverage of oxygen from low to 
high. The leading eigenvalues and corresponding eigenvectors 
at points H and p are plotted in Fig. The inner and outer 
dimensionless radius of the ring are 20 and 30 respectively. 
T = 540K,Po 2 = 1.33 x 10~ 4 mbar; P C o is the bifurcation 
parameter. 




FIG. 4: Bifurcation diagram varying the boundary curvatures 
while fixing ring width. Vi is the local velocity of the prop- 
agating pulse at the inner boundary, vo is the propagating 
speed of a one-dimensional pulse. The increase of bound- 
ary curvature leads to a pulse instability that is associated 
with a Hopf bifurcation, followed again by a turning point. 
T = 540K,Po 2 = 1.33 xl0- 4 mbar,P C o = 4.597 x 10~ 5 mbar. 
r Q — ri is fixed at 8. 



sible for it collapses on the real axis, giving rise to two 
real positive eigenvalues close to point p. Fig. [Sfe) and 
Off) show the corresponding eigenmodes of these two real 
eigenvalues at point p; the absolutely smaller one of the 
two proceeds to cross zero at the turning point. 

For the two-dimensional ring structure, the shape and 
stability of the pulse solutions are dictated by the choice 
of ri and r Q . For more general geometries (e.g. non- 
concentric boundaries) one can parameterize the geome- 
try through the inner curvature, the outer curvature and 
the distance between the boundaries. We will look at 
the effect of geometry taking two distinct one-parameter 
paths: (a) change radii while keeping the ring width 
fixed; and (b) fix the outer boundary and vary the in- 
ner one. 



FIG. 3: The corresponding eigenmodes of three leading eigen- 
values at point H and p as in Fig. |2 (a)-(c) are eigen- 
modes for point H. (a) A = 0.0005, (b) and (c) A = 
—0.0006 db 0.259H. (d)-(f) are eigenmodes for point p. (d) 
A = 0.0026, (e) A = 0.1172, (f) A = 0.2324. 



in the ring, its shape is no longer planar (see the snap- 
shot of a pulse on the stable branch in Fig. [2J; it has to 
curve so as to satisfy the no-flux boundary condition on 
both curved boundaries. The end of the pulse close to 
the inner boundary takes a convex shape and it is known 
that such a convex shape of pulses may lead to instabil- 
ity when the curvature grows above some critical value 
|23q . We will see this boundary effect on the stability of 
pulses more clearly below, when we present bifurcation 
diagrams from a direct continuation in the ring radius. 
After the Hopf bifurcation, the eigenvalue pair respon- 



B. Influence of varying the boundary curvature 

Figure E| reflects the influence of the curvature (of both 
boundaries) on the pulse stability with r Q — n held fixed. 
The leading eigenvalues and corresponding eigenmodes at 
the Hopf bifurcation point are plotted in Fig. [5J Similar 
to the first row of Fig. |3 the eigenmodes corresponding 
to the eigenvalue pair crossing the imaginary axis show 
some structure close to the inner boundary, but here both 
the Hopf and SN bifurcation (at the turning point) are 
induced by the change of boundary curvature alone. The 
Hopf bifurcation is again subcritical (as simulations in- 
dicate). 

In the limit of r^/r = 1 as — > oo, we should have 
an effectively one-dimensional planar pulse with a con- 
stant velocity vq. As we decrease n, the local curvature 
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FIG. 5: Three leading eigenvalues and corresponding eigen- 
modes at the Hopf bifurcation point H in Fig. 0] (a) 
A = -0.0012, (b) and (c) A = -0.0025 ± 0.28991 



of the pulse close to the inner boundary increases and 
the local shape of the pulse becomes more convex (com- 
pare the top inset in Fig. El with a planar pulse). This is 
associated with a decreasing local velocity. For a stable 
two-dimensional pulse expanding in an infinite medium 
as a two-dimensional ring (with azimuthally uniform cur- 
vature), the propagating velocity (and the pulse stabil- 
ity) increases as the ring grows and the pulse becomes 
less convex |22j . There also exists a critical curvature 
above which no pulse solution exists. A pulse "fragment" 
propagating within a two-dimensional ring structure ac- 
quires a non-uniform curvature, so as to satisfy the no- 
flux boundary conditions applied on both curved bound- 
aries (and the linear speed distribution required for it to 
propagate coherently). This can lead to a locally con- 
cave shape of the pulse (larger local velocity) close to the 
outer boundary and a locally convex shape (smaller lo- 
cal velocity) close to the inner boundary. The difference 
in the local curvature of the 2D pulse close to the inert 
boundaries contributes to the spatial structure of the un- 
stable eigenmodes near the Hopf bifurcation point. We 
now study how the pulse responds to changes in V{ while 
keeping r Q fixed. In this case, as we vary r^, both the cur- 
vature of the inner boundary and the distance between 
the two boundaries change. A typical bifurcation dia- 
gram with respect to T{ with fixed r Q is plotted in Fig. El 

When n approaches the outer ring radius r Q , the pulse 
becomes effectively one-dimensional and propagates at a 
limiting velocity of approximately 0.6s -1 (0.02s -1 for uo 
with r Q = 30), consistent with the value read from the 
one-dimensional pulse bifurcation diagram in Fig. ^ The 
qualitative structure of the bifurcation diagram in Fig. [S| 
is similar to that in Fig. 21 

The continuation of pulse solutions with respect to T{ 
for several fixed values of r Q is shown in Fig. {7\ We see 
that, at criticality, the pulse speed at the inner boundary 
(the slowest linear speed along the pulse) does depend on 
the actual value of r Q (and thus the overall pulse shape 
and its nonuniform curvature). All three stable branches 
for different r Q converge to a distinct point at ri/r = 1 , 
that corresponds to the one-dimensional problem (and so 
do, separately, all three unstable ones). 

When the diffusion coefficient is sufficient small, the 
medium can no longer support a pulse solution; this is 
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FIG. 6: Bifurcation diagram of pulses in a 2D ring showing 
the influence of varying the inner ring radius n on the travel- 
ling pulse stability. Point if, marked with a red circle, denotes 
a Hopf bifurcation. Solid (dashed) lines denote stable (unsta- 
ble) steady pulse solutions. For r% close to 30, steady pulse so- 
lutions are computed at discrete points. Insets show the pro- 
files of O coverage. The dimensionless outer radius of the ring 
r Q is fixed at 30. T = 540if , Po 2 = 1.33 x 10 -4 ra6ar, P C o = 
4.597 x lQ- 5 mbar. 
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FIG. 7: Bifurcation diagram varying the inner ring radius n 
while fixing r Q . v\ is the local velocity of the propagating 
pulse at the inner boundary, vo is the propagating speed of 
a one-dimensional pulse. The v\ associated with the Hopf 
bifurcation point depends on r Q . 



associated with a turning point on the continuation curve 
with respect to the diffusion coefficient at fixed n and r Q 
(not shown). Scaling both ring radii by a factor of k 
is equivalent to fixing the geometry but solving a prob- 
lem with a diffusion coefficient scaled by a factor of 1/k 2 . 
This can be used to rationalize the movement of the bifur- 
cation curves towards the right in Fig. Using increas- 
ingly larger r Q should lead to a critical r Q above which 
there exists no pulse solution for a given Vi/r . 
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FIG. 8: Sequences of snapshots showing the decollation of 
propagating pulses from the inner boundary of a ring struc- 
ture, (a) Numerical simulation. The time interval between 
each snapshot is 2 seconds. Posit ion- fixed arrows are used to 
show the change of decollating pulse. The arrows are aligned 
to the boundary of the pulse in the first snapshot. The "tip" 
of the pulse at the inner boundary first slowly moves toward 
the left and shrinkes at the same time, then it decollates from 
the boundary, (b) and (c) are experimental observations. The 
dark area is TiC^ and the gray area is Pt(llO). The white 
propagating pulse is a CO-rich one. Experimental conditions: 
T = 447K, Po 2 = 4 x 10~ 4 ra6ar, (b) P co = 3.35 x 10~ 5 mbar, 
(c) P co = 2.95 x 10~ 5 mbar. 



IV. DYNAMIC SIMULATIONS AND 
EXPERIMENTS 

We now consider how the boundary curvature- 
associated pulse instabilities manifest themselves in time 
for both simulations and experiments. A dynamic sim- 
ulation close to (but just beyond) the Hopf bifurcation 
point shows an initial "tornado like" oscillation of the 
"tip" of the pulse close to inner boundary. This tip be- 
comes increasingly thinner (as indicated by the color con- 
tours) and after some time the pulse "decollates" from 
the inner boundary and eventually dies (Fig. Efa)). 
This transient eventually evolves to the spatially uniform 
steady state which is stable; thus the simulation is consis- 
tent with a subcritical Hopf bifurcation. The limit cycle 
branch is born unstable and "backwards" in parameter 
space; the saddle- type limit cycle and its stable manifold 
constitute the separatrix between the stable pulse and 
the stable uniform solution when these coexist. 

Figure Efb) and c) show this "decollation" mecha- 
nism from experimental observations of a CO-rich pulse 
in a sequence of still shots 0.28 second apart. The di- 
rection of motion of the pulse in the ring is denoted 
by an arrow. The data come from photoemission elec- 
tron microscopy (PEEM) observations of CO oxidation 
on micro-designed Pt catalysts; the black areas are 2000 



A tall Ti layers, deposited through microlithography on 
the Pt(110) single crystal. The inner and outer rings of 
the ring-shaped Pt "corridor" between the tall Ti "walls" 
are 12.5 (160) and 20 (180) /im respectively in Fig. GJb) 
(Fig. Etc)). The experimental conditions are given in the 
figure caption; one must remember, however, when ob- 
serving these experiments that diffusion on the Pt(110) 
crystal is anisotropic. The "thinning" and decollation of 
the pulse close to the inner boundary can be clearly seen. 



V. SUMMARY AND CONCLUSIONS 

We studied the effect of operating conditions and of 
geometry (in particular, of boundary curvature) on the 
speed and stability of pulses in excitable media; our illus- 
trative model described CO oxidation under high vacuum 
conditions on Pt(110). The main instability is a decol- 
lation that occurs at the inner ring boundary; while the 
turning point bounding one-dimensional pulses still exists 
in full two-dimensional studies, the primary destabiliza- 
tion occurs earlier, before the turning point itself. We 
found that this primary instability involves a (subcriti- 
cal) Hopf bifurcation, leading to the decollation of the 
pulse from the inner boundary, and (here) its eventual 
death. 

In microdesigned/microstructured media, pulses often 
encounter sudden changes in curvature (for example, at 
corners, or channel exits); a ring geometry -more pre- 
cisely, families of ring geometries with varying radii- al- 
lows for the systematic study of the effect of steady, 
controlled curvature. This provides insight in the ba- 
sic "boundary decollation" instability which underpins 
many of the events occurring in more complex media. 
Dynamic pattern formation in microchannels is studied 
in many contexts beyond catalytic reactions (see e.g. re- 
cent work in homogeneous reactions |3J , and beyond re- 
acting systems in microti uidics 0, l35j |) . Boundary cur- 
vature and its variations is, we believe, a crucial factor 
in all such applications. Additional issues like boundary 
roughness (statistics of curvature at a much finer scale) 
as well as the use of active boundaries (e.g. consisting 
of a different catalyst, see |3q ) will modify the results 
we described here. We believe, however, that the basic 
reacting front decollation from a curved boundary will 
be a recurring theme in all these contexts, and an impor- 
tant component of pattern formation in complex hetero- 
geneous media. 
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